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Abstract 

A mathematical programming model for a class of single machine family scheduling problems is described 
in this technical report, with the aim of comparing the performance in solving the scheduling problem by means 
of mathematical programming with the performance obtained when using optimal control strategies, that can 
be derived from the application of a dynamic programming-based methodology proposed by the Author. The 
scheduling problem is characterized by the presence of sequence-dependent batch setup and controllable pro¬ 
cessing times; moreover, the generalized due-date model is adopted in the problem. Three mixed-integer linear 
programming (MILP) models are proposed. The best one, from the performance point of view, is a model which 
makes use of two sets of binary variables: the former to define the relative position of jobs and the latter to define 
the exact sequence of jobs. In addition, one of the model exploits a stage-based state space representation which 
can be adopted to define the dynamics of the system. 
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1 Introduction 

A mathematical programming model for the single machine family scheduling problem described in the paper 
“Optimal control strategies for single machine family scheduling with sequence-dependent batch setup and con¬ 
trollable processing times” 11 is here presented in order to compare the dynamic programming-based approach 
proposed in |[T| with a different methodology. Among the possible formulations, a mixed-integer linear program¬ 
ming (MILP) model which makes use of two sets of binary variables (one to define the relative position of jobs and 
the other to define the exact sequence of jobs) is adopted; as a matter of fact, a preliminary performance analysis 
has shown that such a model outperforms other models which solve the same problem. 

This technical report is organized as follows. The scheduling problem is described in section|2]whereas the MILP 
model is reported in section |3] The test of the model (on the numerical example considered in Ul) and an ex¬ 
perimental analysis (on randomly-generated instances of different sizes of the problem) are in section |4] and |5] 
respectively. A comparison of the times required to find an optimal solution between the algorithm proposed in ifTI 
(implemented and solved with Matlab) and the MILP model (implemented and solved with Cplex) is in section |6] 
Two alternative MILP formulations are reported in section |7] (the second of them exploits the stage-based state 
space representation adopted in m). 


2 The scheduling problem 

Consider a single machine where Nk jobs of class Pk, k = 1,..., K, K > 1, have to be executed. All jobs 
belonging to the same class are equivalent. All jobs are available at time instant 0, and preemption is not allowed. 
A sequence of Nk due-dates for jobs of class Pk, namely ddk,i, ddk, 2 , ■ • ■, ddk,Nk^ fc = 1,..., AT, is specihed. It 
is assumed that ddk,i < ddk,i+i for any k = 1,... ,K and i = 1,..., Nk — 1, and that jobs are assigned to the 
due-dates according to the EDD (earliest due-date) rule, in accordance with the generalized due-date model. 

The processing (or service) time of the i-th job of class Pk is a continuous variable ptk,i € k = 

1,... ,K, i = 1,..., Nk, it is assumed that, once the processing time of a job has been chosen, it cannot be 
changed during the service. pt^°’^ is the nominal value of the processing time, which can be compressed up to 
pt^k'^. In accordance with a resource consumption model, the processing time can be expressed as ptk,i{uk,i) = 
p^nom _ Uk,i, with 0 < Uk,i < Uk, Uk = — ptk^)llk, where Uk,i is the amount of a continuous resource 

used to reduce the processing time of the i-th job of class Pk, and jk is the positive compression rate for jobs of 
class Pk- 

A costly setup is required between the execution of two jobs of different classes. In this connection, let stk^h (resp., 
sck,h), k = 1,..., K, h = 1,..., AT, be the setup time (cost) which has to be spent (paid) when switching from a 
job of class Pk to one of class Ph, h f k. No setup is required between the execution of two jobs of the same class 
(i.e., stk^k = 0 and sck,k = 0 V/c = 1,..., AT). Besides, the execution of the first job does not require any setup, 
regardless of the class. Then, stg^k = 0 and sco,fc = 0, k = 1,..., AT, having denoted by sfo,/c and sco,fe the initial 
setup time and the initial setup cost, respectively, when a job of class Pk is the first processed job. 

The performance criterion (cost function) to be minimized is the sum of the total weighted tardiness, the total 
weighted consumption cost (of the continuous resource used to compress the processing time), and the total setup 
cost, namely 

K Nk K Nk K Nk 

EE EE l^k '~ik '^k,i EE ^t^T7k^i,k (1) 

k—1 i—1 k—1 i—1 k—1 i—1 

where: ak,i is the unitary tardiness cost for the i-th job of class Pk', Tk^i is the tardiness of the i-th job of class Pk', 
Pk is the unitary cost related to the deviation from the nominal service time for the execution of jobs of class Pk', 
TTfe i is the index of the class of the job which precedes the Lth job of class Pk in the job sequence, or is the value 
0 if the i-th job of class Pk is the first processed job. The set of values 11 = {tt^ ^ , k = 1,..., K, i = 1,..., Nk} 
is sufficient to define a specific job sequence; moreover, let U = {uk,i , k = 1,..., K, i = 1,..., Nk}. The total 
weighted tardiness depends on both 11 and U, as Tk^i = max{S'fe_i(n, U) -\- i,k — 7fc Uk,i — ddk,i,Q}, 

being S'fe_i(n, U) the time instant at which the Lth job of class Pk starts its execution when values ^ and Uk,i, 
k = 1,..., K, i = 1,..., Nk are adopted. 
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The objective is to find a job sequence If and the values Uk,i, k = ... ,K, i = , Nk, which minimize the 

cost function O- 

3 The mathematical programming model 

3.1 Variables and parameters 

3.1.1 Main decision variables 

• Xh,j,k,i ■ binary variable which defines the relative position of jobs 

- Xhj,k,i = 1 if the j-th job of class Ph is executed before the i-th job of class Pk 

- Xh,j,k,i = 0 Otherwise 

• Sh,j,k,i ■ binary variable which defines the sequence of jobs 

- ^h,j,k,i = 1 if the i-th job of class Pk is the next job after the execution of the j-th job of class Ph 

- ^h,j,k,i = 0 Otherwise 

• Uk,i ■ amount of a continuous used to reduce the processing time of the i-th job of class Pk 

3.1.2 Other decision variables 

• Sk,i ■ time instant at which the f-th job of class Pk starts its execution 

• ptk,i ■ processing time of the f-th job of class Pk 

• T/c i : tardiness of the f-th job of class Pk 

• fife i : setup cost which is paid for the i-th job of class Pk 

• Afe i : setup time which is spent for the z-th job of class Pk 

3.1.3 Parameters 

• Nk : number of jobs of class Pk 

• : nominal value of the processing time of jobs of class Pk 

• : lowest value of the processing time of jobs of class Pk 

• ddk,i ■ due-date for the z-th job of class Pk 

• sth^k '■ setup time which is spent when switching from a job of class Ph to a job of class Pk 

• sch,k '■ setup cost which is paid when switching from a job of class Ph to a job of class Pk 

• Qffc.i : unitary tardiness cost for the z-th job of class Pk 

• j3k ■ unitary processing time’s deviation cost for jobs of class Pk 

• 7 fc : processing time’s compression rate for jobs of class Pk 

• M : large and positive value (“big number” larger than an upper-bound of the maximum completion time) 
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3.2 The MILP model 


3.2.1 Objective function 


EE EE i^k 'Jk ’^k^i H" EE^m 


3.2.2 Constraints 

Tk^i E Sk^i ^k,i “1“ P^kji ddk^i V/u = 1,..., K Vi = 1, . . . , Nk 

K Nh 

^k,i — EE 

h=l j=l 
K Nh 

^k^i — ^ ^ ^ ^ ^ih,k ^hj,k,i ■ • ■ ? ^ — 1 ; • ■ • 5 ^k 

h=l j=l 

K Nh K Nk / K s 


E E E E = E 1 


h=l j=l k=l i=l 
K Nk 


yk = l,...,K Vz=l,...,A^fc 


h=ij=i 
K Nk 


yh = l,...,K = 

fc=li=l 
i 

Y,Sk,^,k,J=0 Vk=l,...,K Vi = l,...,Nk 

i=i 

Nk 

Y.Sk,^,kJ=0 yk = l,...,K \fi=l,...,Nk-2 

j=i+2 


'^k,i E 


ptr“ - 


'ik=l,...,K Vi = l,...,iVfc 


ptk,i = ptr"^ - jkUk,i 'ik = l,...,K yi=l,...,Nk (12) 

Sk,i ^ ^h,j “h ~k pthj ^ (1 ^h,j,k,i) 

Vh = l,...,K yj = l,...,Nh Vk = l,...,K Vz = l,...,A^fe {h,j)^ik,t) (13) 

^h,j ^ ^k,i “I” Ptk,i ^ ^h,j,k,i 

Vh = l,...,K yj = l,...,Nh \/k = l,...,K Vi = l,...,A^fc ih,j)^{k,t) (14) 

Xk,j,k,i = ^ yk = l,...,K Vz = l,...,A^fc Vj = l,...,i-1 (15) 

Xk,j,k,i=0 yk = l,...,K yi = l,...,Nk yj=i,...,Nk (16) 

Xh^j,k,i “1“ — 1 

Vh = l,...,K yj = l,...,Nh \/k = l,...,K Vi = l,...,A^fc {h,j)^{k,i) (17) 


Xh^j,k,i “1“ + Xl ,m,h,j — 2 

'ih = i,...,K 'ij = i,...,Nh 'ik = i,...,K 'ii = i,...,Nk 
{h,j)i^(k,i) (k,i)^(l,m) {l,m)^{h,j) 
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^h,j,k,i ^ 1 ^ (1 






V/i = 1,... 

,K 

Vj = l 


\/k = 

1,.. 

.,K 

Vi = 1 ,.. 

■,Nk (19) 

^h,j,k,i ^ {O 5 1 } 

Vh 

= 1, 

...,K Vj 

= 1,. 

.■,Nh 

Vfc = 1, 

...,K 

Mi 

= 1,. 

..,Nk 

(20) 

^ {0; 1} 

Vh 

= 1,. 

..,K Vj 

= 1,- 

■ .,Nh 

Vfc = 1, 

...,K 

Mi 

= 1,- 

■ .,Nk 

(21) 

Uk,i > 0 

Vfc = 

= 1,.., 

.,K 

Vi = 1,... 

,Nk 







(22) 

Sk,^ > 0 

V/c = 

= 1,.. 

.,K 

Vi = 1,... 

,Nk 







(23) 

ptk,i > 0 

Vfc = 

= 1,.. 

.,K 

Vi = 1,.. 

■ ,Nk 







(24) 

Tk,i > 0 

Vfc = 

= 1,.. 

.,K 

Vi = 1,... 

,Nk 







(25) 

Q.k,i > 0 

Vfc = 

= 1,.. 

.,K 

Vi = 1,.. 

■ ,Nk 







(26) 

> 0 

Vfc = 

= 1,.. 

.,K 

Vi = 1,.. 

■ ,Nk 







(27) 


Constraint (O computes, taking into account (l25T l and the presence of Tk^i in the cost function, the tardiness of the 
i-th job of class Pk which is executed. Constraints (|4]i and (|5]l compute, respectively, the setup cost and the setup 
time which are paid/spent when executing the i-th job of class Pk- Constraint ® ensures that all jobs are executed. 
Constraint (I7]i (resp., ([8]l) guarantees that each job can have one predecessor (successor) only; the inequality is 
necessary to deal with the hrst (the last) executed job. Constraints (|3 and (fTOl i apply the generalized due-date 
model: with respect of class Pk, only jobs i and i + 1 can be (but not necessarily are) consecutive jobs in the 
overall sequence of job executions. Constraint (fTTT i determines, taking into account (l22l l. the interval which is 
allowed for the amount of a continuous resource used to reduce the processing time of the i-th job of class Pk- 
Constraint (fT^ computes the processing time of the of the i-th job of class Pk which is executed. Constraint (fOl l 
(resp., (O) ensures that a job which is scheduled after (before) another job (in accordance with the binary variable 
which dehnes the relative position of jobs) starts actually after (before) that job. Constraints (fTSl l and (fThl l deal 
with the generalized due-date model again: in this case, with respect of class Pk, the constraints ensure that all jobs 
preceding the i-th are scheduled before job i and all jobs following the i-th are scheduled after job r, the redundancy 
of constraints (l9ll- (fT0l i and (flST l- lfThl l in dealing with the generalized due-date model reduces the computational 
times to hnd an optimal solution. Constraints (fTTI) and (fTsT i prevent the system to provide solutions in which some 
jobs constitute a closed/circular subsequence of executions. Constraint (O hxes the values of binary variables 
Xh,j,k,i in accordance with the values of binary variables 6h,j,k,i', as a matter of fact, if the execution of the i-th job 
of class Pk follows the execution of the j-th job of class Ph {Sh,j,k,i = 1) then the j-th job of class Pu dehnitely 
precedes the i-th job of class Pk in the relative position of jobs (xh,j,k,i = 1); such a constraint is equivalent to the 
simpler constraint Xh,j,k,i > Sh,j,k,i, but it provides better performance (lower computational times to solve the 
problem). Finally, constraints (|201 ) - mI271 i dehne the type of the decision variables. 

The proposed problem is a mixed-integer linear mathematical programming problem (MILP). The number of 
binary variables is Nk)^- 


4 Test of the model 

Consider the numerical example for the single machine scheduling problem proposed in IT], in which 4 jobs of class 
Pi and 3 jobs of class P 2 must be executed. The data which characterize the problem are the following. 
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• Due-dates and unitary tardiness costs 


class Pi 

class P 2 

«i,i = 0.75 ddi^i = 19 
0^1,2 — 0.5 ddi^2 — 24 
G!i^3 = 1.5 ddi^^ = 29 
0^1,4 = 0.5 ddi^4 = 41 

cx2,i = 2 dd2^i = 21 
a2,2 = 1 dd2^2 = 24 

a2,3 = 1 dd2,3 = 38 


• Processing time bounds, unitary costs related to the deviation from the nominal processing times, and com¬ 
pression rates 


Setup times and costs 


class Pi 

class P 2 

R 1 pi‘r = 4 , 

/5l 1 p^nom = 8 ^ 

0 1 . pt2" = 4 1 

P 2 1-5 = 6 ^2 1 

)StS 

times 

costs 

sfi,i = 0 Sti ^2 = 1 

SC 14 = 0 sci^2 = 0.5 

St2,l = 0.5 St2,2 = 0 

SC 2,1 = 1 SC2,2 = 0 


By solving such an instance of the MILP problem, the following optimal solution is found. 



1,1 

1,2 

1,3 

{k,i) 

1,4 

1 2,1 

2,2 

2,3 

1,1 

0 

1 

1 

1 

0 

0 

1 

1,2 

0 

0 

1 

1 

0 

0 

1 

1,3 

0 

0 

0 

1 

0 

0 

1 

{h,j) 1,4 

0 

0 

0 

0 

0 

0 

0 

2,1 

1 

1 

1 

1 

0 

1 

1 

2,2 

1 

1 

1 

1 

0 

0 

1 

2,3 

0 

0 

0 

1 

0 

0 

0 



1,1 

1,2 

1,3 

{k,i) 

1,4 

1 2,1 

2,2 

2,3 

1,1 

0 

1 

0 

0 



0 

1,2 

0 

0 

1 

0 



0 

1,3 

0 

0 

0 

0 



1 

{h,j) 1,4 

0 

0 

0 

0 

0 


0 

2,1 

0 

0 

0 

0 

0 

1 

0 

2,2 

1 

0 

0 

0 

0 

0 

0 

2,3 

0 

0 

0 

1 

0 

0 

0 


job {k,i) 

'^k,i 

1,1 

4 

1,2 

0.5 

1,3 

3 

1,4 

0 

2,1 

0 

2,2 

0 

2,3 

0 
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With such values of the main decision variables, the optimal values of the other decision variables are 


job (fc,i) 

Sk,i 

pik,i 

Tk,i 

^k,i 

^k,i 

1,1 

12 

4 

0 

1 

0.5 

1,2 

16.5 

7.5 

0 

0 

0 

1,3 

24 

5 

0 

0 

0 

1,4 

36 

8 

3.5 

1 

0.5 

2,1 

0 

6 

0 

0 

0 

2,2 

6 

6 

0 

0 

0 

2,3 

29 

6 

0 

0.5 

1 


The optimal job sequence provided by the solution is 


J2,l —>• J2,2 Jl,l —>■ Jl.2 Jl,Z —>■ J2,Z Jl,4 

which is the same obtained with the methodology proposed in the paper in the nominal case. Moreover, the total 
cost of the optimal solution is 11.75, which corresponds to the value of the optimal cost-to-go in the initial state 
at to = 0, namely Jq g q(0), provided in |[T1. However, it is worth noting that the optimal values of the decision 
variables (and thus the optimal job sequence) represent an open-loop solution that, when applied, maintains its 
validity as long as no disruption affects the nominal system behavior. Instead, the solution provided in m is a 
closed-loop solution that is able to give the optimal decisions for any actual machine behavior. 


5 Experimental analysis 

An experimental analysis has been carried out to test the performance of the MILP model in solving the considered 
family scheduling problem. Randomly-generated instances of different sizes have been solved. All the parameters 
but (which has been assumed deterministic with 7 fc = 1 Vfc = 1,... ,K) have been considered stochastic 
variables uniformly distributed in the interval [a, b], with a and b reported in the following table. 


parameter 

Uniform distribution [a, b] 
a b 

note 

p^nom 

6 

10 



2 

6 


ddk^i 

ddk.i-i + 0.5 

ddk,i-i -f 12 

ddk.o = 10 

^^h,k 

1 

3 

stk,k = 0 

SCh,k 

0.5 

2.5 

^Ck.^k — 0 

C^k,i 

0.5 

2.5 


13k 

0.5 

2.5 



Only small instances of the problem have been taken into account, which can be solved in seconds or minutes. The 
results of such an analysis follow. 


2 classes - 5 jobs per class 


(Number of problems solved; 100) 


Problem size 


Binary variables 

200 

Other variables 

61 

Constraints 

1227 


Computational time (seconds) 

min 

0.14 

max 

0.86 

mean 

0.29 

stdev 

0.09 
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2 classes - 10 jobs per class 


(Number of problems solved: 100) 


Problem size 


Binary variables 

800 

Other variables 

121 

Constraints 

8757 


Computational time (seconds) 

min 

8.17 

max 

51.39 

mean 

18.91 

stdev 

7.69 


2 classes - 15 jobs per class 


(Number of problems solved: 20) 


Problem size 

Binary variables 

1800 

Other variables 

181 

Constraints 

28587 


Computational time (seconds) 

min 

181.39 

max 

5013.48 

mean 

1466.39 

stdev 

1204.22 


3 classes - 5 jobs per class 


(Number of problems solved: 60) 


Problem size 


Binary variables 

450 

Other variables 

91 

Constraints 

3790 


Computational time (seconds) 

min 

16.36 

max 

232.77 

mean 

42.63 

stdev 

29.03 


Bigger instances of the scheduling problem take several hours to be solved. 


6 Comparison with the dynamic programming-based approach 

The performance of the dynamic programming-based (DP) approach proposed in |[T| and the performance of the 
MILP approach have been compared on some specihc instances of the scheduling problem. 

• DP approach —Implemented and solved with Matlab R2014a 

• MILP approach —Implemented and solved with Cplex 12.6 (with standard settings) 

The results of such a comparison follow. 


2 classes - 5 jobs per class 



Time to find an optimal 
solution (seconds) 




DP Problem size 

DP approach 

0.52 

State space nodes 

61 

MILP approach 

0.25 




MILP Problem size 


Binary variables 

200 

Other variables 

61 

Constraints 

1227 

MILP nodes 

527 
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2 classes - 10 jobs per class 



Time to find an optimal 
solution (seconds) 




DP Problem size 

DP approach 

1.47 

State space nodes 

221 

MILP approach 

11.06 




MILP Problem size 


Binary variables 

800 

Other variables 

121 

Constraints 

8757 

MILP nodes 

13580 


2 classes - 15 jobs per class 



Time to find an optimal 
solution (seconds) 




DP Problem size 

DP approach 

3.17 

State space nodes 

481 

MILP approach 

2903.38 




MILP Problem size 

Binary variables 

1800 

Other variables 

181 

Constraints 

28587 

MILP nodes 

2542617 


2 classes - 20 jobs per class 



Time to find an optimal 
solution (seconds) 




DP Problem size 

DP approach 

5.87 

State space nodes 

841 

MILP approach 

n.a. (> 50000) 




MILP Problem size 


Binary variables 

3200 

Other variables 

241 

Constraints 

66717 

MILP nodes 

n.a. 


3 classes - 5 jobs per class 



Time to find an optimal 
solution (seconds) 




DP Problem size 

DP approach 

4.99 

State space nodes 

541 

MILP approach 

182.84 




3 classes - 10 jobs per class 


MILP Problem size 

Binary variables 

450 

Other variables 

91 

Constraints 

3790 

MILP nodes 

492369 



Time to find an optimal 
solution (seconds) 




DP Problem size 

DP approach 

54.69 

State space nodes 

3631 

MILP approach 

n.a. (> 50000) 




MILP Problem size 


Binary variables 

1800 

Other variables 

181 

Constraints 

28435 

MILP nodes 

n.a. 


4 classes - 5 jobs per class 



Time to find an optimal 
solution (seconds) 




DP Problem size 

DP approach 

91.31 

State space nodes 

4321 

MILP approach 

n.a. (> 50000) 




MILP Problem size 


Binary variables 

800 

Other variables 

121 

Constraints 

8653 

MILP nodes 

n.a. 
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The DP approach outperforms the MILP approach. Only in the case of very small instances of the problem, such as 
2 classes and < 5 jobs, the performance of the two approaches is comparable (both approaches solve the problem 
in less than 1 second), but in all the other cases the DP approach has significantly better performance than the 
MILP approach. 


7 Other mathematical programming models 

Two different mixed-integer linear programming models are proposed in this section. The former is similar to the 
one proposed in subsec tion l3. 21 but it uses only the set of binary variables Sh,j,k,i', the latter exploits the stage-based 
state space representation, as introduced in HI. However, these models do not show better performance (in terms 
of computational times to solve the considered instances of the scheduling problem) than that obtained with the 
model in subsection l3.2l of this technical report. 


7.1 MILP model n.2 

7.1.1 Variables and parameters 

Let, in addition to the variables and parameters introduced in subsection l3.ll 
• Ck,i '■ completion time of the i-th job of class Pk (other decision variable) 


7.1.2 The model 

An alternative MILP model is the following. 

K Nk K Nk K Nk 

min EE EE ^k 'Jk '^k,i “t" (28) 

k—1 i—1 k—1 i—1 k—1 i—1 

subject to 

Ck,i = Sk,i + Ak,i + ptk,i \/k = l,...,K \/i = l,...,Nk (29) 


Tk,^ > Ck,i - ddk,^ yk = l,...,K yi = l,...,Nk (30) 


K Nh 

^ ^ ^ ^ V/l = 1, . . . , K Vi = 1, . 

■ .,Nk 

(31) 

h=l i=l 

K Nh 

Ak^i — ^ ^ ^ ^ dhj^k,i V/c = 1, . . . , K Vi = 1, . 

...Nk 

(32) 

h=l3=1 

K Nh K Nh . K . 

E E E E = (E 1 


(33) 

h=l3 = lk=li=l ^k=l ' 

K Nh 

Y.Y.^h,3,k,i<l Vfc = l,...,iT 'il=l,...,Nk 


(34) 

h=l3 = l 

K Nk 

EE ^ 1 V/i = i,...,Ar 'ij = i,...,Nh 


(35) 

k—li—1 

i 

^4,i.fcj=0 'ik = l,...,K Vi = l,...,iVfe 


(36) 


i=i 
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Nk 


Y,Sk,^,k,J=0 yk = l,...,K \/i=l,...,Nk-2 

(37) 

j=i+2 


uk,^<— - yk = i,...,K yi = i,...,Nk 

Ik 

(38) 

ptk,i = - lkUk,i 'ik = l,...,K 'ii = l,...,Nk 

(39) 

Sk,i > Ch,j — M (1 — Sh,j,k,i) 



Vh = l,...,K yj = l,...,Nh Vk = l,...,K yi = l,...,Nk {h,j)^{k,i) (40) 

K Nh 

Vfc=l,...,X V^ = l,...,iVfc (41) 

h=lj=l 

W,ze{o,i} yh = i,...,K yj = i,...,Nh yk = i,...,K yi = i,...,Nk (42) 


Uk,i > 0 

Vfc = 1,. 

..,K 

Vi = 1,. 

■ ■,Nk 

(43) 

Sk,^ > 0 

Vfc = 1,. 

..,K 

Vi = 1,. 

...Nk 

(44) 

Ck,^ > 0 

Vfc = 1,. 

..,K 

Vi = 1,. 

...Nk 

(45) 

ptk,i > 0 

Vfc = 1,. 

..,K 

Vi = 1,. 

.■,Nk 

(46) 

Tk,i > 0 

Vfc = 1,. 

..,K 

Vi = 1,. 

■ ■,Nk 

(47) 

rikj > 0 

Vfc = 1,. 

..,K 

Vi = 1,. 

.■,Nk 

(48) 

> 0 

Vfc = 1,. 

..,K 

Vi = 1,. 

.■,Nk 

(49) 


Constraint (l29T l computes the completion time of the i-th job of class Pk which is executed. Constraints (I30l l-f- (l39l l 
have the same meanings of constraints ©-MO described in subsection 13.21 Constraint (l40l i ensures that the *-th 
job of class Pk is scheduled after the j-th job of class P^ when 6h,j,k,i = 1- Constraint (HTt determines the comple¬ 
tion time of the hrst job which is executed (since the hrst job is the only job for which it is ^^=1 ^f=i = 

0). Finally, constraints (|421 i-m|49]| dehne the type of the decision variables. 

The proposed problem is a mixed-integer linear mathematical programming problem (MILP). The number of 
binary variables is 


7.2 MILP model n.3 

The scheduling model can be represented through a state space model in which the system state, when a new 
decision has to be taken, namely at time instant tj, j = 0,1,..., N — 1, with N = is the {N + 2)- 

dimensional vector Xj = [nij,..., nK,j, being Ukj, k = 1,... ,K, the number of jobs of class Pk 

already completed at time instant tj, and hj the class of the last served job. As a matter of fact, decision instants 
are discrete in time; they correspond to the initial time instant and the instants at which a job is completed, but for 
the last one. 

It is apparent that the system state does not change between two successive decision instants. The {N + l)-tuple 
{nij,..., riKj, hj) will be referred to as the discrete part of the system state vector; it can be represented through 
a stage-based state transition graph in which the states at the j-th stage, j G {0,1,..., At — 1}, are those such that 
SfeLi = j- An example of a state transition graph is reported in Fig. [T] for the case of two classes of jobs, 
with A^i = 4 and N 2 = 3. 

The system state evolves with the application of the control = f{Xj) at the decision instants fj, j = 0,1,... ,N— 
1. The control, which is a function of the system state, is relevant to the choice of the class of the next job to be 
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Figure 1: State transition graph of the discrete part of the system state (in the case K = 2, with = 4, and 
N2 = 3) 


served and the value of the processing time. Then, Wj = [(5i j ,..., Skj, where 6k,j £ {0,1}, fc = 1,..., K, 
is a binary decision variable whose value is 1 if a job of class Pk is selected for the (j + l)-th service, and 0 other¬ 
wise, and Tj is the processing time for the selected job. Obviously, it must be X]fc=i ^k,j = IV j = 0,1,.. .,7V — 1. 
Moreover, it it worth noting that the choice of considering the processing time Tj as a decision variable instead 
of the amount of the continuous resource used to compress the processing time, namely Uj, does not modify the 

matter of the problem, as Uj = — - — 7 — ( Y^k=i P^k°^ ^k,j — Tj). 


The state equations for the considered model are 



ni,j+i 


J + ^1,1 

2^3 + 1 = 

nK,j+i 

= 

nK,j + 6 k, j 


Tij+i 


^k=lj^ ^k,j 


. ^3 + 1 . 


- ^3 + Tj + J2k=l ^'^hj,k 6k,j . 


(50) 


and the initial state vector is Xq = [0,..., 0,0, 0]"^. 

The objective is to find values Sk,j, k = 1,... ,K, and Tj, j = 0,. .., N — 1, that 

Af-l 

max{fj+i — + Phj+i {p^hj+i ~ 

3=0 


minimize the objective function 


Tj) + SChj,hj+i 


(51) 


7.2.1 Variables and parameters 

Let, in addition to the variables and parameters introduced in section lLTl 
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• Xk,i,j : binary variable whose value is 1 if the i-th job of class Pk is executed at the j-th stage, and 0 otherwise 

(main decision variable). It is apparent that j (then variable can be employed in place 

• Tj : processing time of the job which is executed at the j-th stage (main decision variable) 

• Vlj : setup cost which is paid at the j-th stage (other decision variable) 

• Kj : setup time which is spent at the j-th stage (other decision variable) 

• Sj : time instant at which the job which is executed at the j-th stage starts its execution (other decision 
variable) 

• Cj : completion time of the job which is executed at the j-th stage (other decision variable) 


• Ck,i '■ completion time of the r-th job of class Pk (other decision variable) 

• N : total number of jobs (parameter) 

7.2.2 The model 

An alternative MILP model is the following. 

K Nk K Nk N-1 

min ^ ^ ak,i Tk^i + X! X! “ Ptk,i) + ^ % (52) 

k—1i—1 k—1i—1 0 

subject to 

Tk,^>Ck,^-ddk,^ 'ik = l,...,K Vi = l,...,iVfc (53) 

ptk,i > pt'k'^ Vfc = 1,..., AT Vi = 1,..., (54) 

ptk,^<ptT'" yk = l,...,K yi = l,...,Nk (55) 

> sch,k i^Xh,ij-i+^Xk,ij - Ij yj = l,...,N-l yh,k = l,...,K (56) 

tlo = 0 (57) 

/ tVh Nk , 

> sth,k i^Xh,t,j-i+^Xk, 2 ,j - l] 'ij = 1,...,N -1 yh,k = l,...,K (58) 

^ i=l ^ 

Ao = 0 (59) 

Vj = l,...,iV-l (60) 

^0 = 0 (61) 

Cj = Sj +Aj+Tj Vj = 0,..., iV - 1 (62) 

Tj > ptk,i - M {1 - Xk,^,J) Vj =0,...,iV-l \fk=l,...,K Vi = l,...,iVfc (63) 

Sk,^ > Sj - M {1 - Xk,^,J) Vj = 0,...,iV-l Vfc = l,...,iV yi = l,...,Nk (64) 

Ck,^ > Cj - M (1 - Xk,i,j) yj = 0,...,N-l yk = l,...,K Vi = l,...,iVfc (65) 

Sk,^>Ck,^-l yk=l,...,K Vi = 2,...,iVfc (66) 
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K Nk 

=1 Vj = 0,...,iV-l (67) 

k—1i—1 
N-l Nk 

Vfc=l,...,X (68) 

N-l 

Y^Xk,^,,=l Vk=l,...,K Vi = l,...,iVfc (69) 

j=o 

xk,^,Je{0,l} yk = i,...,K v* = i,...,iVfe yj = o,...,N-i (70) 

r, >0 Vj=0,...,iV-l (71) 

Sk,^>0 yk = l,...,K yi = l,...,Nk (72) 

Ck,^>0 Vfc = l,...,X yi = l,...,Nk (73) 

ptk,^>o yk = i,...,K yi = i,...,Nk (74) 

Tk,,>0 yk = l,...,K yi = l,...,Nk (75) 

ilj > 0 Vj = 0,..., - 1 (76) 

A, >0 Vj = 0,...,lV-l (77) 

Sj>0 Vj = 0,..., iV - 1 (78) 

Cj>0 \/j = 0,...,N -1 (79) 


Constraint (l53T l computes, taking into account (fTsT l and the presence of ^ in the cost function, the tardiness of the 
i-th job of class Pk which is executed. Constraints (l54l i and (fSST l set the range of allowed values of the processing 
time of the i-th job of class Pk. Constraints (l56]) and (fSST l compute, respectively, the setup cost and the setup 
time which are paid/spent at the j-th stage. Constraints dSTT i and (l59l l set the initial setup cost and the initial setup 
time, respectively, to 0. Since it can be shown that an optimal solution of the considered single machine family 
scheduling problem is characterized by the absence of idle times between the execution of two subsequent jobs 
(even if the cost function is not of regular type), constraint (l60l) states that each job starts as soon as the previous 
one completes, and constraint dMl l sets the start time of the first job to 0. Constraint (l62l) computes the completion 
time of the job which is executed at the j-th stage. Constraint (|63| sets the processing time of the job which is 
executed at the j-th stage, taking into account the values of binary variables Xk,i,j as well as the range fixed by (l5^ 
and (1551) . Constraints (l64l) and (1651) determine, respectively, the start time and the completion time of the i-th job of 
class Pk which is executed. Constraint (l66]) deals with the generalized due-date model as it ensures that, for what 
concerns class Pk, the i-th job is scheduled after the {i — l)-th job (that is, the r-th job starts after the completion 
time of the {i — l)-th job). Constraint (l67T i guarantees that a job is executed only once. Constraint dhST l ensures that 
all jobs are executed, whereas constraint (l69l l establishes that, in each stage, only one job can be served. Finally, 
constraints (l70l) ^(l79ll define the type of the decision variables. 

The proposed problem is a mixed-integer linear mathematical programming problem (MILP). The number of 
binary variables is {N,k=i ^k)^- 
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